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ABSTRACT 

Disk accretion to rotating stars with complex magnetic fields is investigated using 
full three-dimensional magnctohydrodynamic (MHD) simulations. The studied mag- 
netic configurations include superpositions of misaligned dipole and quadrupolc fields 
and off-centre dipoles. The simulations show that when the quadrupolc component 
is comparable to the dipole component, the magnetic field has a complex structure 
with three major magnetic poles on the surface of the star and three sets of loops of 
field lines connecting them. A significant amount of matter flows to the quadrupolc 
"belt" , forming a ring-like hot spot on the star. If the maximum strength of the mag- 
netic field on the star is fixed, then we observe that the mass accretion rate, the 
torque on the star, and the area covered by hot spots are several times smaller in 
the quadrupolc-dominant cases than in the pure dipole cases. The influence of the 
quadrupolc component on the shape of the hot spots becomes noticeable when the 
ratio of the quadrupolc and dipole field strengths B q /B<i > 0.5. It becomes dominant 
in determining the shape of the hot spots when B q /Bd > 1. We conclude that if the 
quadrupolc component is larger than the dipole one, then the shape of the hot spots 
is determined by the quadrupole field component. In the case of an off-centre dipole 
field, most of the matter flows through a one-armed accretion stream, forming a large 
hot spot on the surface, with a second much smaller secondary spot. The light curves 
may have simple, sinusoidal shapes, thus mimicking stars with pure dipole fields. Or, 
they may be complex and unusual. In some cases the light curves may be indicators 
of a complex field, in particular if the inclination angle is known independently. We 
also note that in the case of complex fields, magnctosphcric gaps are often not empty, 
and this may be important for the survival of close-in exosolar planets. 

Key words: accretion, accretion disks - magnetic fields - MHD - stars: magnetic 
fields. 



1 INTRODUCTION 

The magnetic field of a rotating star can have a strong in- 
fluence on the matter in an accretion disk. The field can 
disrupt the disk and channel the accreting matter to the 
star along the field lines. The associated magnetic activ- 
ity can be observed through photometric and spectral mea- 
surements of different types of stars, like young solar-type 
Classical T Tauri stars (CTTSs) (Hartmann et al. 1994), X- 
ray pulsars and millisecond pulsars (Ghosh & Lamb 1978; 
Chakrabarty et al. 2003) , cataclysmic variables (Wickramas- 
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inghe et al. 1991; Warner 1995, 2000), and also brown dwarfs 
(e.g., Scholz & Ray 2006). The topology of the stellar mag- 
netic field plays an important role in disk accretion and in 
the photometric and spectral appearance of the stars. 

The dipole magnetic model has been studied since 
Ghosh & Lamb (1979a, b), both theoretically and with 2D 
and 3D MHD simulations. However, the intrinsic field of 
the star may be more complex than a dipole field. Safier 
(1998) argued that the magnetic field of CTTSs may be 
strongly non-dipolar. Zeeman measurements of the magnetic 
field of CTTSs based on photospheric spectral lines show 
that the strong (1 — 3 kG) magnetic field is probably not 
ordered, and indicates that the field is non-dipolar close to 
the star (Johns-Krull et al. 1999, 2007). Other magnetic field 
measurements with the Zeeman-Doppler imaging technique 
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have shown that in a number of rapidly rotating low-mass 
stars, the magnetic field has a complicated multipolar topol- 
ogy close to the star (Donati & Cameron f997; Donati et 
al. 1999; Jardine et al. 2002). Recently, observations from 
the ESPaDOnS/NAVAL spectropolarimeter (Donati et al. 
2007a) revealed that the magnetic field geometry of CTTSs 
BP Tau, V2129 Oph and SU Aur is more complex than 
dipole. 

Theoretical and numerical research has been done on 
accretion to a rotating star with a non-dipole field. Jar- 
dine et al. (2006) investigated the possible paths of the 
accreting matter in the case of a multipolar field derived 
from observations. Gregory et al. (2006) developed a sim- 
plified stationary model for such accretion, von Rckowski 
and Brandenburg (2006) did axisymmetric simulations, and 
von Rekowski and Piskunov (2006) did 3D simulations of 
the disk-magnetosphere interaction in the case where the 
magnetic field is generated by dynamo processes in both the 
star and the disk. They obtained a time-variable magnetic 
field of the star with a complicated multipolar configuration 
which shows that the dynamo may give rise to a complex 
magnetic field structure. Recent observations and analysis of 
the magnetic field and brightness distribution on the surface 
of CTTSs suggest that the magnetic field is predominantly 
octupolar with strength B oct ~ 1.2 kG and a much smaller 
dipole component B d ~ 0.35 kG (Donati et al. 2007b). Most 
of the accretion luminosity is concentrated in a large spot 
close to the magnetic pole. Such observations are interest- 
ing and can be compared with different 3D MHD simulation 
models. Here we investigate, using 3D MHD simulations, 
disk accretion to rotating stars with dipole plus quadrupole 
fields and off-centre dipole fields. 

In our earlier work we performed full 3D simulations 
of disk accretion to stars with pure quadrupole and aligned 
dipole plus quadrupole fields (Long, et al. 2007) and found 
that the funnel streams and associated hot spots have dif- 
ferent features compared with the case of accretion to a star 
with pure dipole field. When the quadrupole component is 
significant, matter always accretes to a quadrupole "belt" 
and forms a ring-like hot spot on the surface of the star. 

In this paper we show the results of full 3D MHD sim- 
ulations of accretion to stars with more complex magnetic 
fields. Compared with Long et al. (2007), we consider the 
more general case where the dipole and quadrupole are mis- 
aligned (see §3). Also, we consider the case of accretion to 
a star with internal point dipoles displaced from the star's 
center (see §3.4). Further, we performed a set of simulations 
with the aim of understanding the properties of stars with 
different contributions from the dipole and quadrupole com- 
ponents. 

One point of interest is to determine the how the area 
of the star's surface covered by hot spots depends on the 
the quadrupole/dipole ratio. A further point of interest is to 
determine how the torque on the star depends on the struc- 
ture of the field. That is, we can compare the torque for 
cases with a complex field with those having a dipole field. 
We discuss all these points in §4. We calculate the light 
curves and discuss the use of the light curves to investigate 
the structure of the magnetic field of the star. The conclu- 
sions and discussion of possible applications of our results 
are given in §5. Some high resolution figures and animations 
are available at http://www.astro.cornell.edu/~long. 



2 THE NUMERICAL MODEL AND THE 
MAGNETIC CONFIGURATIONS 

2.1 Model 

In a series of previous papers (Koldoba et al. 2002; Ro- 
manova et al. 2004; Ustyugova et al. 2006, Long, et al. 2005, 
2007), we described the model used in our 3D MHD simu- 
lations. The model used in this paper is similar to that of 
Long et al. (2007) and therefore it is only briefly summarized 
here. 

We consider a rotating magnetized star surrounded 
by an accretion disk with a low-density, high-temperature 
corona above and below the disk. The disk and the corona 
are initially in a quasi-equilibrium state. We solve the 3D 
MHD equations in a reference frame rotating with the star, 
with the z— axis aligned with the star's rotation axis. The 
magnetic field is decomposed into a "main" component B 
and a variable component B± . Here, B is time-independent 
in the rotating reference frame and consists of the dipole and 
quadrupole parts: Bo = Bd + B q . Further, i?i is in general 
time dependent due to currents in the simulation region and 
is calculated from the MHD equations in the rotating refer- 
ence frame. 

A special "cubed" sphere grid with the advantages of 
both spherical and cartesian coordinate systems was devel- 
oped (Koldoba et al. 2002). It consists of N r concentric 
spheres, where each sphere represents an inflated cube. Each 
sphere consists of 6 sectors corresponding to the 6 sides of 
the cube, and an N x N grid of curvilinear cartesian coordi- 
nates is introduced in each sector. Thus, the whole simula- 
tion region consists of six blocks with N r x N 2 cells. In the 
current simulations we chose a grid resolution of 75 x 31 2 . 
Other resolutions were also investigated for comparison. The 
coarser grids give satisfactory results for a pure dipole field. 
However, for a quadrupole field, the code requires a finer grid 
because of the higher magnetic field gradients. The code is 
a second order Godunov-type numerical scheme (see, e.g., 
Toro 1999) developed earlier in our group (Koldoba et al. 
2002). Viscosity is incorporated into the MHD equations in 
the interior of the disk so as to control the rate of matter in- 
flow to the star. For the viscosity we used an a— prescription 
with a — 0.04 in all simulation runs, and with smaller/larger 
values for testing. 

Initial Conditions: The region considered consists of the 
star located in the center of coordinate system, a dense disk 
located in the equatorial plane and a low-density corona 
which occupies the rest of the simulation region. Initially, 
the disk and corona are in rotational hydrodynamic equilib- 
rium. That is, the sum of the gravitational, centrifugal, and 
pressure gradient forces is zero at each point of the simu- 
lation region. The initial magnetic field is a combination of 
dipole and quadrupole field components which are force-free 
at t = 0. The initial rotational velocity in the disk is close 
to, but not exactly, Keplerian (i.e., the pressure gradient is 
taken into account). The corona at different cylindrical radii 
r rotates with angular velocities corresponding to the Ke- 
plerian velocity of the disk at this distance r. This initial 
rotation is assumed so as to avoid a strong initial disconti- 
nuity of the magnetic field at the boundary between the disk 
and corona. The distribution of density and pressure in the 
disk and corona and the complete description of these initial 
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conditions is given in Romanova et al. (2002) and Ustyugova 
et al. (2006). 

The initial accretion disk extends inward to an inner 
radius r d and has a temperature T d which is much less than 
the corona temperature Td = 0.01T C . The density of the 
corona is 100 times less than the density of the disk, p c — 
O.Olpd. These values of T c , T d , p c , p d are specified at the disk- 
corona boundary near the inner radius of the disk. 

Boundary Conditions: At the inner boundary (r = R — 
R*), where R* is the radius of the star, boundary conditions 
are applied to the density dp/dr — 0, pressure dp/dr — 
0, entropy dS/dr = 0, velocity d(v - fi x R)/3r = 
and the magnetic field, dBe/dr — 0, dB^/dr = 0. The 
r— component of the magnetic field satisfies d(r 2 B r )/dr = 0. 
The matter flow is frozen to the strong magnetic field so that 
we have (v — $7 x R) || B. At the outer boundary, free bound- 
ary conditions are taken for all variables with the additional 
condition that matter is not permitted to flow in through 
the outer boundary. The investigated numerical region is 
large, ~ 45-R*, so that the initial reservoir of matter in the 
disk is large and sufficient for the performed simulations. 
Matter flows slowly inward from the external regions of the 
disk. During the simulation times studied here only a small 
fraction of the total disk matter accretes to the star. 
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Table 1. Sample reference units for typical CTTSs, cataclysmic 
variables, and millisecond pulsars. Real dimensional values for 
variables can be obtained by multiplying the dimcnsionlcss values 
of variables shown in this paper by these reference units. See also 
§2.2. 



2.2 Reference Units: 

We solve the MHD equations using dimensionless vari- 
ables: distance R = R/Ro, velocity v = v/vo, time t — 
t/Po, etc. The subscript "0" denotes a set of reference (di- 
mensional) values for variables, which are chosen as fol- 
lows: _Ro = R*/0.35, where 7?* is the radius of the star; 
v = (GM/R ) 1/2 ; time scale P = 2ttR /v . Other ref- 
erence values are: angular velocity Qq = vo/Ro; magnetic 
field B = -B»o(-R*/i?o) 3 , where B»o is the reference mag- 
netic field on the surface of the star; dipole magnetic mo- 
ment po = BqRq; quadrupole moment Do = BoRq; den- 
sity po = Bq/vo; pressure po = po«o> mass accretion rate 
Mo = poVoRo; angular momentum flux Lo = pov 2 R%\ en- 
ergy per unit time Eg = poVoR 2 (the radiation flux J is also 
in units of Eo, see §3.1); temperature To = IZpo/po, where 
1Z is the gas constant; and the effective blackbody tempera- 
ture T e ff,o = (poVo/c) > where a is the Stefan-Boltzmann 
constant. 

In the subsequent sections and figures, we show dimen- 
sionless values for all quantities and drop the tildes (~). To 
obtain the real dimensional values of variables, one needs to 
multiply the dimensionless values by the corresponding ref- 
erence units. Our dimensionless simulations are applicable 
to different astrophysical objects with different scales. For 
convenience, we list the reference values for typical CTTSs, 
cataclysmic variables, and millisecond pulsars in Tab. 1. 



2.3 Magnetic Configurations 

Combination of Dipole and Quadrupole Fields: The 

intrinsic magnetic field of the star is B = —Vtp, where the 
scalar potential of the magnetic field is <p(r) = Em a /|r — r a |, 
m a is the magnetic "charge" analogous to the electric 
charge, and r and r a are the positions of the observer and the 
magnetic "charges" respectively. The scalar potential can be 



represented as a multipole expansion in powers of 1/r and 
the quadrupole term is 

(p (2> = D al3 n a n/3/2r :i , 

where n a — x a /r, x a are the components of r, D a p is the 
magnetic quadrupole moment tensor, summation over re- 
peated indices is implied, and D aa = 0. In the axisymmetric 
case where flu = D22 = — D33/2, let D — D33 be the value 
of the quadrupole moment and refer to the axis of symme- 
try as the "direction" of the quadrupole moment, so we have 
the quadrupole moment D. The combination of dipole and 
quadrupole magnetic fields can be written as: 



B (r) = 3 ^f^ + ^ ( 5 ( D.^-l )r - 



3D 

2T 4 



(E>. ?)£>,(!) 



where r and D are the unit vectors for the position and 
the quadrupole moment respectively. In general, the dipole 
and quadrupole moments pL and D are misaligned relative 
to the rotational axis CI, at angles O and 9d respectively. 
In addition, they can also be in different meridional planes 
with an angle $ between the fl — pi and fl — D planes. 

The values of p and D determine the magnetic field 
strength on the star's surface. For example, if /x and D are 
aligned with the rotational axis ft, and p, = 0.5 , D = 0.5, 
then the strength of the dipole component on the north pole 
of the star is B,d = 0.68 kG, and the quadrupole one is 
B, q = 0.32 kG. 

Off-centre Magnetic Configurations: The magnetic 
field of a star may deviate from the center of the star. Such 
displacement naturally follows from the dynamo mechanism 
of magnetic field generation. If the dipole and quadrupole 
moments are displaced by a and b from the center of the 
star, we can replace r with (r — a) and (r — b) in the 
corresponding terms in Eqn. 1 to obtain the magnetic field 
around the star. In this paper we investigate accretion to a 
displaced dipole and to a set of displaced dipoles. 
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Figure 1. The distribution of dipole and quadrupole components 
of magnetic fields. The left-hand panel represents the case when 
fi = D = 0.5, © = 0°,©o = 45°; the right-hand panel represents 
the case when fi = D = 0.5, 6 = 30°, © D = 45°, * = 90°. The 
strengths of the quadrupole and dipole are comparable at r ~ 1 
where the disk is stopped by the magnetosphere. 

3 ACCRETION TO STARS WITH DIFFERENT 
COMPLEX MAGNETIC FIELDS 
CONFIGURATIONS 

In Long et al. (2007), we have shown results of accretion to 
stars with dipole plus quadrupole fields with aligned axes. 
In this section, we present simulation results of accretion 
to a star with misaligned dipole plus quadrupole magnetic 
fields. First, we consider the configuration when the dipole 
magnetic moment \i is aligned with the stellar spin axis fi, 
that is, = 0°, but the quadrupole magnetic moment D 
is inclined relative to the spin axis at an angle ©_d = 45°. 
Next, we consider a more general configuration when /x, D 
and fi are all misaligned relative to each other, and there 
is an angle $ between the (fi — /i) and (fi — D) planes. 
To ensure that the case is general enough to represent a 
possible situation in nature, we choose O = 30°, 9d = 45° 
and $ = 90°. 

The strengths of the dipole and quadrupole moments 
are the same in both cases, jx = D — 0.5. Here we choose a 
quite strong quadrupole component to get an example of a 
complex magnetic field, so that the disk is disrupted not by 
the dipole component only but by the superposition of the 
dipole and quadrupole fields. Fig. 1 shows the distribution 
of magnetic fields for these two cases in +x direction. We 
can see that the quadrupole component is still comparable 
to the dipole component at r ~ 1 where the disk is stopped 
by the magnetosphere (see Fig. 5). This means that both 
dipole and quadrupole contribute to stop the disk. 

3.1 Misaligned dipole plus quadrupole 
configurations 

Now we consider the configuration when the (fi — fi) and 
(fi - D) planes coincide: /j, = D = 0.5, = 0°, Q D = 45°. 
Fig. 2 shows the strength of the magnetic field on the sur- 
face of the star in different projections. One can see that the 
field is not symmetric and shows three strong magnetic poles 
and a much weaker magnetic pole with different polarities on 
the surface, due to the asymmetry of the misaligned dipole 
plus quadrupole configuration. The left-hand panel shows a 
strong positive magnetic pole (red) with B = 67 (dimen- 
sionless value) and a strong negative pole (dark blue) with 
B — —41. The middle panel shows a second strong posi- 
tive pole (yellow) with B — 35, and in fact a very weak 
and small negative pole in the middle of the surface with 
B — — 7 (dark green). We should note that the strength 




Figure 2. The strength of the magnetic field on the surface of the 
star for the case fi = 0.5, D = 0.5, = 0°, © D = 45° at different 
projections: from the equatorial plane (left-hand panel), close to 
the yz plane (middle panel), and from the south pole (right-hand 
panel). The red and blue regions represent different strength and 
polarities. The strength varies from the most blue point B = —41 
to the most red point B = 67 (dimcnsionless value, see §2.2 and 
Tab. 1) and reaches B = at the boundary between the yellow 
and green areas. 




Figure 3. The complex magnetic field lines for the case fi = 0.5, 
D = 0.5, 6 = 0°, e D = 45°, at t = (top panel) and t = 8 
(bottom panel). The strength of the magnetic field along the field 
lines varies from red (positive magnetic pole) to blue (negative 
pole). 



of the magnetic field around this weak pole is more nega- 
tive than at the pole. The right-hand panel shows both the 
strong negative pole (dark blue) which is extended and the 
second strong positive pole (yellow). 

Fig. 3 shows a three-dimensional (3D) view of the struc- 
ture of the magnetic field lines. The color shows the differ- 
ent strengths and polarities of the magnetic field lines, from 
maximum positive (red) to maximum negative (blue) val- 
ues. The magnetic field lines are more complex than in the 
aligned dipole plus quadrupole case. Now we see three sets 
of loops of closed field lines in the region near the star, con- 
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Figure 4. 3D views of matter flow to the star for the case p = 0.5, 
D = 0.5, = 0°, Br, = 45°, at t = 8. The top panel shows the 
view from the north, the middle panel from the south. The green 
surface is a constant density surface, p = 0.27 (dimensionless 
value, see §2.2). The bottom panel shows the equatorial plane. 
The colors of the magnetic field lines correspond to those in Fig. 
3. Only the inner part of the simulation region is shown. 



necting the three major poles and the weaker pole on the 
star, which is different from the one big and one small loops 
of field lines shown in the aligned dipole plus quadrupole 
case (Long et al. 2007). 

Fig. 4 shows a 3D view of matter flow around the star. 
The disk matter is disrupted by the complex field and is 
lifted above the equatorial plane at the magnetospheric ra- 
dius r = r m , where the magnetic stress balances the matter 
stress, /3 — (p + pv 2 )/(B 2 /8tt) = 1. We can see that most 
of the matter flows between the loops of field lines to the 
extended negative magnetic pole by choosing the shortest 
path which is energetically favorable, forming a modified 
quadrupole "belt". This "belt" has a different shape from 
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Figure 6. The hot spots viewed at different angles: from the 
rotational equatorial plane (left-hand panel); from near the yz 
plane (middle panel); and from the south pole (right-hand panel). 
The red color corresponds to the densest region with maximum 
density p = 1.4 (dimensionless value, see §2.2). The bottom panel 
shows the light curves for different inclination angles. 

that in the case of aligned dipole plus quadrupole config- 
urations (see Long et al. 2007), which looks like a more 
"regular" sheet perpendicular to the magnetic axis. Here, 
the "belt" is twisted, although it is approximately perpen- 
dicular to the quadrupole axis. The bottom panel shows a 
slice in the equatorial plane, and we can see that in the —x 
direction, the matter penetrates between field lines to the 
surface of the star. 

Fig. 5 shows different projections of the accretion flow. 
One can see that the matter flow is not symmetric in the 
xz and xy planes. The xz projection shows that in — x di- 
rection, the magnetic stress is weak due to the combination 
of the dipole and tilted quadrupole components, so that the 
matter is not lifted, and flows almost directly to the sur- 
face of the star, and the equatorial plane does not have the 
magnetospheric gap which is typical for pure dipole cases. 
Because the quadrupole moment is inclined at 45° relative 
to the z— axis in xz plane, in the +x direction, there is 
one set of closed loops of magnetic field lines, which stops 
the disk and forms a magnetospheric gap. Comparison of 
the dipole and quadrupole components in the +x direction 
shows that in the region where the disk stops (/3 = 1), the 
dipole and quadrupole components are approximately equal, 
so that the quadrupole has a strong influence on the matter 
flow, which is clearly seen from the above figures. In the yz 
projection, the matter flow is relatively symmetric because 
SI, fi and D are all in the xz plane. In the xy projection, we 
can see that there is a gap in +x direction where there is no 
direct matter flow, but there is no gap in the —x direction. 

Fig. 6 shows the hot spots on the surface of the star 
at t = 8. One can see that there is a ring-like hot spot and 
a small round hot spot nearby. There is no clear hot spot 
near the positive part of the D and fi axes, due to the 
weak funnel stream in the northern hemisphere. The ring- 
like hot spot comes from accretion to the quadrupole "belt" , 
which is approximately symmetric relative to D because the 
quadrupole magnetic component B q is stronger than the 
dipole component Bd in the region close to the star. A small 
round spot forms as a result of the connection of field lines 
between the weak magnetic pole (B — —7) and the nearby 
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Figure 5. Disk accretion viewed from different projections for the case of fi = 0.5, D = 0.5, © = 0°, ©£> = 45°. The left-hand panel 
shows the projection in xz plane, the middle panel in the yz plane, and the right-hand panel in the xy plane. The color background 
shows the density distribution varying from p = 0.01 (blue) to p = 1.9 (red), and the magnetic field lines are shown in red. The yellow 
lines correspond to j3 = 1. 



region where B is positive. Part of the matter flows to the 
south and forms this round hot spot. 

We calculated the light curves from the hot spots on 
the rotating star's surface. Assuming the total energy of the 
inflowing matter is radiated isotropically as blackbody ra- 
diation, we obtain the flux of the radiation in a direction 
k, 

J = r 2 F obs = J I(R, k) cos 8dS, (2) 

where r is the distance between the star and the observer, 
Fobs is the observed flux, I(R, k) is the specific intensity of 
the radiation from a position R on the star's surface into a 
solid angle element dfi in the direction k, 9 = arccos (R ■ k), 
AS is an element of the surface area. The specific intensity 
can be obtained as 

I(R,k) = -F e (R)cos6, (3) 

TV 

where F e (R) is the total energy flux of the inflowing matter. 
In our simulations, J is the received energy per unit time, 
and the dimensionless value of J is in units of Eq which 
is shown in Tab. 1 and discussed in §2.2. The hot spots 
constantly change their shape and location. However, our 3D 
simulations show that the changes are relatively small. So 
we choose the hot spots at some moment of time, fix them, 
and rotate the star to obtain the light curves. The bottom 
panel of Fig. 6 shows the light curves at t — 8 for different 
inclination angles i = arccos (O ■ k). One can see that the 
light curve is approximately sinusoidal for small inclination 
angles. This is because the observer can only see a part of 
the ring-like hot spots at this time and the sinusoidal shape 
is determined by the rotation of the spots with the star. One 
can see that the shapes of the light curves at i = 60° and 
90° are unusual and do not correspond to any pure dipole 
field case (Romanova et al. 2004). When i increases, the 
small round hot spot and more of the ring-like hot spots 
can be observed, and consequently the peak intensity and 
variability of the light curves are larger. 

3.2 A more general case: fj,,D,fl all misaligned 

Next we investigate the more general case when the dipole 
moment fi, quadrupole moment D and spin axis ft of the 
star are all misaligned relative to each other: fi = 0.5, D = 



0.5, 9 = 45°, D = 30° and $ = 90°. Fig. 7 shows the 
distribution of the magnetic field on the surface of the star. 
We can again see that there are three dominant magnetic 
poles with magnetic strengths B = 65, B = 37andB — —44 
respectively. There is another weak magnetic pole in the 
middle region shown in the middle panel, but it shrinks and 
partially merges into the boundary of the nearby pole and 
looks smaller than that in the previous case. Fig. 8 shows 3D 
plots of matter flow at t = 8. In the top panel, one can see 
that again, a significant amount of matter flows through the 
quadrupole "belt" below the loops of field lines, because the 
quadrupole component dominates near the star. And a single 
strong stream flows to the region nearby the north pole. So 
it could be expected that ring-like hot spots and round hot 
spots form on the surface of the star. The bottom panel 
shows a slice in the xy plane. One can see that some matter 
can flow to the star directly without leaving the equatorial 
plane. We also can see three sets of loops of field lines. 

Fig. 9 shows different projections of the matter flow. 
The matter flow is no longer symmetric in any projection 
which is expected due to the magnetic configuration. In the 
xy and yz planes, we also can see that there are three sets of 
loops of field lines which regulate the shape of the accretion 
flow. In the northern hemisphere, some matter is lifted from 
the equatorial plane and flows to the high-latitude area near 
the dipole axis. This is because the accretion disk is close 
to one of the major magnetic poles in this direction. In the 
southern hemisphere, matter penetrates the /3 line between 
loops of field lines, and forms a modified quadrupole "belt". 
The magnetosphere is not empty in the xy plane, and matter 
can go directly to the star in some directions. 

Fig. 10 shows the hot spots and the associated light 
curves. We can see there is one strong arc-like hot spot near 
the dipole magnetic axis, but not at the magnetic pole. This 
hot spot comes from the one-armed stream shown in Fig. 9. 
Another ring-like hot spot is below the equatorial plane but 
not symmetric with respect to either the fl, fi or D axis. 
The light curve is sinusoidal for the small inclination angle 
i = 15°. At larger i beside the round north spot, part of the 
ring-like spot becomes visible to the observer, so the light 
curves depart from the sinusoidal shapes and have a larger 
amplitude. 
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Figure 9. Slices of density distribution for a more general case (p = 0.5, D = 0.5, = 30°, &]j = 45°, <J> = 90°). The left-hand, middle 
and right-hand panels show the projection in the xz, yz and xy planes respectively. The color background shows the density distribution 
varying from p = 0.01 (blue) to p = 2.1 (red). The magnetic field lines are shown in red. The yellow lines shows where /3 = 1. 




Figure 7. Distribution of the magnetic field strength on the star's 
surface for the case p = 0.5, D = 0.5, = 45°, D = 30° and 
$ = 90° at different projections: from the equatorial plane (left- 
hand panel), from near the yz plane (middle panel), and from the 
south pole (right-hand panel). Red and blue colors show positive 
and negative polarities. The magnetic strength varies from B = 
-44 to B = 65. 



3.3 Off-centre dipole fields 

Until now, all we discussed are cases in which the dipole 
and quadrupole magnetic moments are located at the cen- 
ter of the star. The interior convective circulations may be 
displaced from the center and lead to off-centre magnetic 
configurations. Therefore we now consider disk accretion to 
a star with an off-centre pure dipole field, with /j, = 0.5. The 
magnetic moment is located at x — 0.5R,, y — 0, z — 0, 
that is, it is shifted from the centre of the star by half ra- 
dius of the star. The misalignment angle is Q = 30°. Fig. 11 
shows the 3D plot of disk accretion for this case. One can see 
that the north and south magnetic poles are not symmetric 
on the surface of the star. Due to the shift and the tilt of 
the dipole moment, the north magnetic pole is closer to the 
accretion disk than the south pole, so that most of the mat- 
ter flows to the north magnetic pole. Fig. 12 shows the hot 
spots on the surface of the star and the corresponding light 
curves. The light curves are sinusoidal at i = 15°, 30°, and 
have unusual shapes at i = 60° , 90° . 

In another set of runs, we chose a superposition of sev- 
eral dipoles placed at different places in the star, and with 
different orientations of their axes. Fig. 13 illustrates the 
considered configurations of three off-centre dipoles of equal 
value fi = 0.5, displaced by 0.4i?„ from the center of the star 
in the equatorial plane and azimuthally separated by 120°. 
In configuration (a), all moments are misaligned relative to 
the rotation axis at O = 30°, and in configuration (b), they 
are arranged in the equatorial plane in an anti-clockwise 
manner. 




Figure 8. 3D views of disk accretion to a star in the more general 
misaligned dipole plus quadrupole case at t = 8, when p = 0.5, 
D = 0.5, = 45°, O d = 30° and $ = 90°. A constant density 
level (p = 0.25, dimensionless value) of the disk is shown in green 
in the top panel; different density levels in the equatorial plane are 
shown in the bottom panel. The field lines are shown in multiple 
colors. 



Fig. 14 shows the distribution of the magnetic field on 
the surface of the star. One can see that there are six poles 
on the star. In case (a), three positive and three negative 
poles are in northern and southern hemispheres respectively. 
In case (b), all poles are in the equatorial planes and poles 
with different polarities alternate. 

Fig. 15 shows a 3D view of the matter flow for the cases 
(a) and (b). In case (a), most of the matter flows to the 
north poles which are closer to the accretion disk due to 
the inclination of the dipoles. The magnetospheric gap is 
empty in the equatorial plane. In case (b), all the matter 
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Figure 10. Hot spots from different angles at t = 8 for the case 
of p = 0.5, D = 0.5, 6 = 45°, e D = 30° and $ = 90°. The 
left-hand, middle and right-hand panels represent edge-on, top 
and bottom views respectively. The red color corresponds to the 
densest region with maximum density p = 2.3. The bottom panel 
shows the light curves for different inclination angles. 

flows between loops of the closed magnetic field lines in the 
equatorial plane and forms some interesting equatorial fun- 
nels. 

Fig. 16 shows the hot spots for the cases (a) and (b). 
There is a big triangular hot spot around the north pole for 
the case a, which spans the area near the magnetic poles. 
Another hot spot is present near the south pole but with 
very low densities. This confirms the asymmetric matter flow 
shown in Fig. 15. We can also see that there are no polar hot 
spots for case (b). Only ring-like hot spots are present in the 
equatorial plane, which corresponds to what we observed in 
Fig. 15. 



4 CONTRIBUTION OF THE QUADRUPOLE 
TO PROPERTIES OF MAGNETIZED STARS 

In this section we analyze properties of magnetized stars at 
different ratios between the dipole and quadrupole fields. 
First, we consider different configurations at the same max- 
imum value of the field on the surface of the star. Then 
we analyze different properties such as the area covered by 
hot spots, mass accretion rate and spin torque. Next, we 
fix the dipole component but change the contribution of the 
quadrupole component with the main goal of understanding 
how strong the quadrupole should be compared with the 
dipole in order to change the shape of the hot spots from 
the pure dipole cases. 

4.1 Area covered by hot spots, and torque 

In this section we compare properties of magnetized stars of 
different configurations under the condition that the max- 
imum magnetic field on the surface of the star is the the 
same in all cases. 

We choose a number of configurations ranging from 
purely dipole to purely quadrupole: (a) pure dipole field, 
with /x = 1.07, = 45°; (b) pure quadrupole field, with 
D = 0.5, 6d = 0°; (c) aligned dipole plus quadrupole field, 
with /i = 0.5, D = 0.27, 9 = Od = 0°; (d) misaligned 




Figure 11. 3D views of matter flow to a star with an off-centre 
dipole at t = 8, where p = 0.5, = 30° and the dipole moment 
is located at x = 0.5i?*,y = 0, z = 0. The green color in the top 
panel shows one of the density levels p = 0.3. The bottom panel 
shows the equatorial plane with different density levels. The field 
lines are shown in multiple colors. 




Figure 12. Top panel: hot spots for the off-centre case shown in 
Fig. 11, at t = 8. The red color corresponds to the densest region 
with maximum density p = 2.6, and the yellow color in the right- 
hand panel represents the density p = 0.7. Bottom panel: the 
light curves. 



dipole plus quadrupole field, with fi = D = 0.37, O = 0°, 
On = 45°. In all cases, the maximum strength of the surface 
magnetic field is the same, about B max = 50 (dimensionless 
value), but the contributions of the dipole and quadrupole 
components are different. In cases (a), (b) and (c), the field 
has a maximum value at the north pole, and in case (d), it 
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Figure 15. 3D views of disk accretion to a star with the three off-centered dipoles described in Fig. 13, at t 
(a); Right-hand panel: case (b). 



Left-hand panel: case 





Figure 13. The configurations of three off-centre, equally sepa- 
rated dipole magnetic moments, which are 0.4iJ» away from the 
center of the star in the equatorial plane, with (J-i = P2 = = 
0.5. (a) The magnetic moments are tilted away from the rotation 
axis at = 30°; (b) the magnetic moments are arranged in the 
equatorial plane in a circular fashion. 



is near the quadrupole axis. Here, for case (a), we choose a 
relatively large O to avoid MHD instabilities (see Kulkarni 
& Romanova 2008, Romanova, Kulkarni & Lovelace 2008). 
The discussed configurations are shown in Fig. 17. Below, 
we compare different properties of accretion to these stars. 



o o 



Figure 14. The surface magnetic field on a star with the three 
off-centre dipoles described in Fig. 13. Top panel: case (a); bottom 
panel: case (b). 



field decreases faster with distance than the dipole field, so 
the the disk comes closer. 



4-1.1 Where does the disk stop? 

We calculated the magnetospheric radius r m in +x direction 
in the equatorial plane and obtained the value of r m for 
above four cases: r m , a = 1.47, r m>b = 0.62, r m . c = 1.0, 
r m ,d = 1-05 respectively. We can see that the disk stops 
closer to the star for the pure quadrupole case than for the 
pure dipole case, and r m has intermediate values for the 
mixed cases. This result is expected because the quadrupole 



4-1.2 Area covered by hot spots 

We also calculated the area covered by hot spots in the above 
cases (a)-(d). Like in Romanova et al. (2004), we chose some 
density p and calculated the area A(p) of the region where 
the density is larger than p. Then we obtain the fraction of 
the stellar surface covered by the spots, f(p) — A(p)/A*. 

We also consider the temperature distribution on the 
surface of the star. In spite of the possible complex processes 
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Figure 18. Left-hand panel: area covered by the hot spots where 
the density is larger than p (in units of po, see §2.2); right-hand 
panel: the area covered by the hot spots where the effective tem- 
perature is larger than T e g (in units of T c ff q, see §2.2). The solid, 
dashed, dotted and dash-dotted lines represent the four cases (a)- 
(d) shown in Fig. 17 respectively. 



Figure 16. The hot spots from different angles for the cases de- 
scribed in Fig. 13, at t = 8. Top panel: case (a), the maximum 
densities in the hot spots in the northern and southern hemi- 
spheres are p = 1.4 (red) and p = 0.1 (light blue) respectively. 
Bottom panel: case (b), the maximum density is p = 0.8. 



(a) 




(b) 



(c) 



(d) 




Figure 17. The configurations of different magnetic fields: case 
(a), pure dipole field; case (b), pure quadrupole field; case (c), 
aligned dipole plus quadrupole field; case (d), misaligned dipole 
plus quadrupole field. 



of radiation from the hot spots, we suggest that the total 
energy of the incoming stream is radiated approximately 
as a blackbody. The total energy flux carried by inflowing 
matter to the point R on the star is, 

F e (K) = pn-v(^v 2 +w), (4) 

where n = — f , v is the matter velocity, and w = r y(p/ p) /("■/— 
1) is the specific enthalpy of the matter. So we have F e (R) = 
<tT% and the effective blackbody temperature, 



on ■ v , 1 o n 
a 1 



1/4 



(5) 



where a is the Stefan-Boltzmann constant. Similarly, we ob- 
tain the area A(T c a) with effective temperature larger than 
T e a, and the fraction f(T eS ) = A(T eB )/A*. 

Fig. 18 shows the distributions of f(p) and /(T e ff) for 
cases (a)-(d). We can see that the area covered by spots in 
the pure quadrupole case (case b) is several times smaller 
than that in the dipole case (case a). The area covered by 
spots in the mixed dipole plus quadrupole cases (b and c) 
is in between. We conclude that if most of the field is de- 
termined by the quadrupole, the expected area covered by 
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Figure 19. The mass accretion rate M (left-hand panel) and 
the angular momentum fluxes associated with the magnetic fields 
Lj (right-hand panel) for stars with different magnetic configura- 
tions. The solid, dashed, dotted and dash-dotted lines represent 
the four cases (a)-(d) shown in Fig. 17 respectively. 



spots is smaller. We do not know exactly the reason for 
the smaller fraction / in the cases of quadrupole field. One 
of the factors which we should mention is that the mass 
accretion rate to the surface of the star is also smaller in 
the case of quadrupole configurations (see §4.1.3 and Fig. 
19). It seems that it is "harder" for matter to penetrate 
between the quadrupolar field lines to the poles compared 
with the dipole case, and matter accretes through a narrow 
quadrupole "belt" (Long et al. 2007), which leads to smaller 
accretion rates and smaller hot spots. Fig. 18 also shows 
that for equal f(p) or f(T e g), the density and temperature 
of the spots in the pure quadrupole case are smaller than in 
the pure dipole case. 

4-1.3 Mass accretion rate and torque 

As we mentioned above, matter accretion rate to the surface 
of the star M = — J dS ■ pv p is several times smaller for the 
quadrupole case than for the dipole case (v p is the poloidal 
component of the velocity). The mixed cases (c) and (d) 
have intermediate M (see Fig. 19). We also calculated the 
angular momentum fluxes transported from the accretion 
disk and corona to stars for cases (a)-(d). The angular mo- 
mentum flux carried by the matter L m — — J dS ■ prv$v p is 
about 10 times smaller than the angular momentum flux as- 
sociated with the field Lf = J dS ■ rB^Bp/ (4tt). So we only 
show Lf in Fig. 19. When matter flows in, the magnetic field 
lines inflate from the initial configuration. Later, they recon- 
nect, and the transport of angular momentum flux becomes 
strong. This process is longer in the pure dipole case (a). For 
corotation radius r cor = 2 considered in the above cases, the 
star spins up mostly through the field lines connecting it to 
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the inner part of the disk. Here we choose slowly rotating 
stars (r cor — 2) for all magnetic configurations to make sure 
that the magnetospheric radius r m is smaller than the coro- 
tation radius r cor , and the star spins up in all cases. So we 
can compare the positive spin torques for different magnetic 
geometries. Fig. 19 shows M and Lf for cases (a)-(d). One 
can see that the angular momentum flux is largest in the 
pure dipole case (a) and smallest in the pure quadrupole 
case (b). This is an interesting result because in the case of 
the quadrupole field, the disk comes closer to the star and 
it rotates faster, so that one would expect larger angular 
momentum flux in the case of the quadrupole field. But we 
see the opposite. We can speculate that stars with predomi- 
nantly multipolar fields may not spin up as fast as stars with 
a mainly dipole field. This can be explained by the smaller 
connectivity between closed multipolar field lines and the 
inner region of the disk (see also Long et al. 2007). 



4.2 Hot spot shape 

Next, we fix the dipole magnetic field at some relatively 
high value /i = 2, the misalignment angle at some value 
= 30°, and gradually increase the quadrupole component, 
D — 0.5, 1, 1.5, 2. In all cases, the quadrupole is aligned with 
the fl axis. The main goal is to understand the ratio D/n 
at which the properties of hot spots will depart from the 
dipole ones, and the ratio D/[i at which the quadrupole will 
strongly influence the shape of hot spots. 

Fig. 20 shows that in the pure dipole case (left column), 
the hot spots have a typical arc-like shape. With the in- 
crease of the quadrupole component, the southern hot spot 
is stretched. When the quadrupole is strong enough, the 
southern hot spot forms a ring. We conclude that the influ- 
ence of the quadrupole component on the shape of the hot 
spots becomes noticeable when the ratio of the quadrupole 
and dipole moments D/fj, > 0.25, and becomes dominant in 
determining the hot spot shape when D/fi > 0.5. The corre- 
sponding ratios of the magnetic field strengths on the surface 
of the star are B q /Bd = 0.5 and B q /Bd = 1 respectively. 

Fig. 21 shows the light curves for the above cases. It 
is interesting to note that at i = 30° the amplitude is rel- 
atively small for the strong quadrupole case. This can be 
explained by the fact that part of the southern ring-like hot 
spot is not visible to the observers in this orientation. When 
the quadrupole component become more important, the in- 
tensity and amplitude of the light curves become smaller. 
The shape of the light curves is more irregular than in the 
pure dipole configuration. Light curves are complex but may 
provide information about the magnetic configuration of the 
star, if the inclination angle is known independently. 



5 SUMMARY AND DISCUSSION 

We investigated disk accretion to rotating stars with dif- 
ferent complex magnetic fields using full three-dimensional 
magnetohydrodynamic simulations. The main results of this 
work are the following: 

1. We investigated accretion to a star with a dipole plus 
quadrupole field for general conditions where the moments 
of the dipole and quadrupole components are misaligned 
and also where they are in different meridional planes. We 




Figure 20. The hot spots for different cases at t = 8. The top row 
shows one side of the star, and the bottom row the other. Each 
column represents different configurations D = 0, D = 0.5, D = 
1.0, D = 1.5 and D = 2.0 from the left to the right respectively. 
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Figure 21. Light curves for inclination angle i = 30° (top panel) 
and i = 60° (bottom panel) at t = 8. The solid, dashed, dash- 
dotted, dotted and dash-dot-dotted lines represent the configura- 
tions D = 0, D = 0.5, D = 1.0, D = 1.5 and D = 2.0 respectively. 



concentrated on the case where both dipole and quadrupole 
fields are strong so that the total field is complex. In this 
case, there are three major magnetic poles with different 
polarities on the surface of the star, and three sets of loops 
of closed magnetic field lines connecting these poles. The 
matter flow is not symmetric; it tends to flow to the star 
along the shortest path to the nearest magnetic pole. Most 
of the matter flows to the star through a narrow quadrupole 
"belt" between loops of field lines. Thus arc-like and ring- 
like hot spots typically form on the star. 

2. We investigated cases with first one and then a set 
of off-centre dipoles. In the case of one displaced dipole, the 
matter flow is not symmetric. More matter flows through one 
stream than the others, and one hot spot is larger than the 
others. In the case of three displaced dipoles, there are three 
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positive and three negative poles on the star. Depending on 
the directions of the dipole moments, matter flows through 
multiple streams and forms strong hot spots at the north 
pole or in the equatorial plane. 

3. The mass accretion rate, the area covered by hot 
spots, and the torque on the star are several times smaller 
in the quadrupole-dominated case than in the dipole cases 
(for conditions where the maximum strength of the magnetic 
field on the star's surface is the same). 

4. The quadrupole component has a noticeable influ- 
ence on the shape of hot spots if the magnetic field of the 
quadrupole component B q > 0.5Bd {D/n > 0.25). The 
shape of hot spots is chiefly determined by the quadrupole 
field if B q > B d (D/fi > 0.5). 

5. The light curves from hot spots of rotating stars with 
complex magnetic fields are often sinusoidal in the case of 
small inclination angles, i ^ 30°, but are more complex for 
i ^ 60°. Even for i ^ 60°, the light curves often resemble 
those of misaligned dipoles at large i. However, in some cases 
they are very unusual. We conclude that, (a) very unusual, 
non-sinusoidal light curves may be sign of the complex fields; 
(b) simple sinusoidal light curves do not rule out complex 
fields; (c) light curves may be a tool for analyzing the com- 
plexity of the field, if the inclination angle is determined 
independently. 

Many CTTSs show highly variable light curves. From 
our simulations we conclude that hot spots do not change 
their shape or position significantly and thus may not be 
responsible for such strong variability. The variability might 
be explained by some other mechanism, such as a highly 
variable accretion rate, or unstable accretion (Romanova, 
Kulkarni & Lovelace 2008; Kulkarni & Romanova 2008). 

Results of 3D MHD simulations of accretion to stars 
with complex fields can be used to compare models with ob- 
servations. For example, recently, Donati et al. (2007b) de- 
rived the magnetic topology on the surface of CTTS V2129 
Oph, and found that it has a dominant octupole compo- 
nent with B oct ~ 1.2kG and a weaker dipole component 
with Bd ~ 0.35kG. Most of the accretion luminosity is con- 
centrated in a quite large high-latitude spot (5% of the to- 
tal stellar surface) close to the magnetic pole. We have not 
modelled accretion to stars with octupole fields, but with 
different contributions from dipole and quadrupole compo- 
nents. This gives information on how matter flows to the 
star. Our analysis shows that the quadrupole dominantly 
determines the matter flow close to the star and the shape 
of hot spots if B q > Bd- Thus, we would expect that in a 
star with a dominant octupole, the accretion spot would be 
determined by one of octupole belts, or part of the belt in 
the case of tilted octupole fields. So from our point of view, 
for such a strong octupole field, the single round spot in the 
high latitude region is an unexpected feature, which requires 
additional analysis in the future. 

The results of our simulations are also important for 
analysis of magnetospheric gaps in young stars. In stars 
with multiple magnetic poles, for example, in the case of 
the dipole plus quadrupole fields with different orientations 
of the axes, the probability is high that one accretion stream 
crosses the equatorial plane, thus filling it with matter. This 
factor may be important in determining the survival of close- 
in exosolar planets, which have a peak in their spatial dis- 
tribution at a few stellar radii. In the case of a dipole field, 



a large low-density magnetospheric gap may be formed in 
the equatorial plane, which may halt subsequent migration 
of planets (Lin et al. 1996, Romanova & Lovelace 2006). 
However, if a star has a complex magnetic field with several 
poles, one of streams may flow through the equatorial plane 
and the magnetospheric gap may be not empty. 
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